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1. Introduction 

Classical Yang-Mills fields have been the subject of intense study in the literature. 
Besides being interesting objects in their own right, much of the motivation arises from 
the relevance of non-abelian gauge theories for our understanding of the interactions 
among elementary particles. Although, in that case, one is concerned with quantum 
fields, particular classical configurations have been argued to play a fundamental role in 
explaining some phenomena. In this context, the pioneer work of Polyakov [l[ showed 
how pseudoparticle solutions become relevant in the semiclassical approximation to 
the path integral. These solutions are then argued to be responsible of some of the 
most intriguing effects of Quantum Field Theory. This is the case for the Confinement 
property of the 3-dimensional compact abelian model, explained originally by Polyakov. 
In 4 dimensions 't Hooft showed that instantons || provide an explanation for the 
so-called [/(l)-problem. This triggered a joint effort of physicists and mathematicians 
in the 70's, which culminated with the construction of all self-dual solutions to the 
euclidean Yang-Mills equations compactified on a sphere [|J. On physical grounds, 
the compactification on a sphere is a condition equivalent to the requirement of finite 
classical action. In other words, one is looking for solutions which exist isolated, i.e. 
surrounded by the perturbative vacuum (F^ u = 0). However, in Quantum Field Theory 
neither the requirement of finite action nor the one of isolated solutions is compulsory. 
For example, it was suggested by the work of Saviddy and the Copenhagen group 
that the QCD vacuum was permeated by some structures carrying chromo-magnetic 
flux. These are neither isolated nor carry finite action. Henceforth, it could also be 
relevant and interesting to study classical solutions with different boundary conditions. 



In this respect a good deal of interest has focused on the study of Yang-Mills fields 
on the 4- dimensional torus (for a review see Ref. 0). The torus adds topological 
features which have an appealing physical interpretation ||. Furthermore, solutions 
on the torus can be considered solutions on R 4 which are periodic, with a total action 
which diverges (with finite action per unit cell). Unfortunately, despite the effort, only 
very special self-dual solutions on the torus, having constant field strength, are known 



analytically || [K| . The lattice formulation of Yang-Mills theory has proven a precise 



and efficient method for obtaining these solutions numerically jO], . A basic building 
block of many of these configurations is a certain lump carrying fractional topological 
charge Q = 4. These lumps cannot exist isolated — their size is determined by the 
distance to their neighbours. The moduli space of a self-dual solution on the torus, 
given by the index theorem, has four parameters per lump, which can be associated to 
the coordinates of its center. This suggests that one can obtain self-dual solutions on 
R 4 by deforming these configurations away from the periodic arrangement (at no cost 
of action). A picture of the Yang-Mills quantum vacuum as a liquid of these fractional 
topological charge lumps, as suggested by our group |T3), could account at the same 
time for the observed string tension and topological susceptibility. 

On the analytical side, an interesting tool is provided by the Nahm transforma- 
tion flnfl . This transformation maps a U(N) self-dual gauge field on the torus with 
topological charge Q into a U (Q) self-dual gauge field on the dual torus with topolog- 
ical charge N. The use of this transformation, which can be considered a particular 
case of a duality transformation, allows one to prove that, in the absence of twist, 
there are no Q = 1 instanton solutions on the torus (see for example ||15|| ). Its use has 
proved useful as well for constructing new instanton solutions in S\ x R 3 flTjJ. For 
the torus case, it has been shown that the constant field strength solutions in SU(2) 
are mapped into themselves JT^]. In this paper we will show that it is possible to use 
lattice techniques to explicitly construct the Nahm transform of a given self-dual lat- 
tice gauge field configuration. This opens the door to the possibility of systematically 
investigating the properties of the Nahm transformation for the whole set of self-dual 
gauge field configurations. In the paper we will apply the method to the case of one of 
the non-trivial self-dual solution on the torus, producing accurate results for its trans- 
form. As a matter of fact, this together with the study of the zero modes for the Dirac 
equation in the adjoint representation, which are supersymmetric partners of the gauge 
field, provides new numerical descriptions of the torus self-dual gauge fields, with some 
relative advantages over the conventional discretization. There is also hope that, as 
in the aforementioned S\ x R 3 case, this study could eventually lead to an analytical 
approach to these classical fields. 

The layout of the paper is as follows. In the next section, we present the details 
of the Nahm transformation and the main formulas which relate our problem to the 
solution of homogeneous and inhomogeneous Dirac equations. In the following section, 
we explain how one can formulate the latter on the lattice, and describe the numerical 



techniques that we use to solve them. In section 4, we apply the technique to the actual 
construction of the Nahm transform for some input self-dual configurations. We first 
test the method by looking at its result in a known case, and then apply it to other 
non-trivial self-dual configuration. Finally, in section 5 we summarize our results and 
discuss future prospects. 



2. The Nahm transformation 

Let us consider a 4 dimensional torus of size lo x l± x l 2 x l 3 , and let / M represent the 
vector (0, ... ,1^, ... ,0), whose only non-zero component is the ji th component. Now 
consider a self-dual gauge field configuration A^x) defined on this torus. It satisfies: 

A v (x + ij = %(x)A v (x)n+(x) + int l (x)d v n+(x) , (2.1) 

where Q^x) are the twist matrices. For SU(N), these matrices must fulfill the consis- 
tency condition: 

Q^x + v)Q v {x) = Z^ u Q„(x + iijtt^x) . (2.2) 

When fields transforming in the fundamental representation of the gauge group appear, 
all the constants must be equal to 1; otherwise, they are in general elements of 
the center of the group Z^, which can be parametrized as = exp(2mn^ u /N). The 
twist tensor n^ u is antisymmetric, and its elements are integers defined modulo N. 

Let Q stand for the topological charge of the gauge field configuration. As is well- 
known, the Atiyah-Singer index theorem implies that the difference between the number 
of positive chirality and negative chirality solutions of the Dirac equation for fermion 
fields transforming in the fundamental representation of the gauge group is given by Q. 
Let us assume that for our gauge configuration there are no negative chirality solutions. 
Then, there are exactly Q positive chirality solutions, which we will label ^f a (x), with 
a = 1, . . . , Q. They satisfy: 

BV a (x) = , (2.3) 

where D = D^T M is the positive chirality Weyl operator and = <9 M — iA^. In the 
Weyl basis we have: 

where T M = (/, —la) and F M = r+. Furthermore, the solutions satisfy the following 
boundary conditions: 

* a {x + g = %(x)V*(x) . (2.5) 
Now consider the family of gauge fields: 

A^x,z) = A 11 {x) + 2ttz^I, (2.6) 



where z^ are 4 real numbers. For all z, the field strength is the same, and hence 
they are all self-dual and have the same topological charge. Therefore, we obtain a 
family ty a (x,z) of positive chirality solutions of the Dirac equation: 

D z ty a (x,z) = (D -2mzl)y a (x,z) = , (2.7) 
satisfying the boundary condition Eq. |275| and normalized as follows: 

d 4 x (^ a (x,z)) + yP(x,z) = S af3 . (2.8) 



Now notice that exp(— 2mz fi x ft ) ^ a (x, z + z) satisfies the same equation than ty a (x, z). 
However, in general, the boundary conditions are different, since the right-hand side of 
Eq. ^75] gets multiplied by exp(—2mz fi l fi ). This new factor becomes simply unity if z^ 

is an integer multiple of 1/7 M . Hence, defining the vector Z M = (0, . . . , j-, . . . , 0) we can 
write: 

V a (x, z + t) = tf^z, z){n'+{z))„ a exp(2mxM . (2.9) 

This is so because any solution can be written as a linear combination of the basis 
functions tyP(x,z). The coefficients (z)) p a cannot in general be chosen equal to 1, 
if we insist in i&P(x, z) being continuous in z. 

Now let us construct the Nahm transform of the gauge field A^(x). It is given by: 

(A fl (z)) a/3 = i J d 4 x (* a (x,z)) + £-* P (x,z) . (2.10) 

This is a U(Q) gauge field defined on the dual torus (of size ^ x ~ x x ~). Using 
Eq. one finds that A^(z) satisfies a relation analogous to Eq. ^TT] (exchanging the 
roles of x and z) in terms of Q' . Now, one can in terms of this field construct the field 
strength tensor F^ u . The Nahm-transformed gauge field has the following properties: 

• F^v is again self-dual. 

• The first and second Chern classes and the ranks for the original and transformed 
gauge fields are related through: 

rk{F) = c 2 {F) - l -c\{F) (2.11) 



Cl (F) = - / {dzp A dx^Y A c x {F) (2.12) 

c 2 (F)=rk(F) + ~cl(F) . (2.13) 

Thus, in particular, one sees that the roles of the rank of the group and the 
topological charge are exchanged by the Nahm transformation when the first 
Chern class vanishes. 
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From the previous statement it follows that if we start with SU(N) gauge fields 
(with no twist = mod N), then the Nahm transform is in SU(Q). 



The Nahm transformation is an involution: if we apply it twice we go back to the 
original gauge field. Thus, it can be considered a duality transformation. 

If we start with a family of gauge fields depending on some parameters, then the 
Nahm transformation generates a new set of self-dual gauge fields depending on 
those parameters. Hence, we have induced a mapping between the moduli spaces 
of the gauge field and its transform. This mapping is an isometry with respect 
to the natural metric of these moduli spaces. 



For a proof of these properties, see Ref. [15 



Now, let us consider the vicinity of a point z in the dual torus. We can make a 
Taylor expansion of the positive chirality solutions of the Dirac equation in the vicinity 
of this point: 

V a (x, z + Az) = V a (x, z) + Az^(x, z) + Az^Az^ix, z) + . . . (2.14) 

By plugging this equation into the Dirac equation for z + Az and equating powers of 
Az^ on both sides, we obtain for the first two orders: 

D z ^ a (x,z) = (2.15) 
D z %(x,z) = 2mT^ a (x,z) (2.16) 

Now, in terms of these functions, and defining: 

Pf{z) = = J d 4 x(^(x,z)) + ^(x,z) 

Q$(z) = <*g*ft (2.17) 

one can write the vector potential and the field-strength tensor coming out of the 
Nahm transformation as follows: 

= iP^z) (2.18) 
F»v{z) = % (Q^(z) - Q+ V (z) + [P^z), P v {z)\) (2.19) 

The normalization conditions imply that A^z) is hermitian. Henceforth, to obtain the 
Nahm-transformed gauge field at one point z of the dual torus, one has simply to solve 
Eqs. |2.15H2.16| . None of the two equations has a unique solution. Choosing a solution 



within both sets of equations amounts to a choice of gauge for the Nahm-transformed 
fields. More specifically, any normalized set of solutions of Eq. |2.15| is related by a 
unitary transformation Q' a Jz) to any other. Furthermore, once this choice is made, it 



can be easily shown that selecting a particular solution \[f"(x,z) of Eq. |2.16j amounts 



to the choice of a gauge in the neighbourhood of z. One can obtain any solution of 



Eq. p,16| by adding to general solution of Eq. |2.15| : ^(x, z) + S^ a z), 



where must be antihermitian due to the normalization conditions. This produces 
a change in the vector potential — > A^ + iSft, but as can be readily verified from 
the equations above, F^ u is left invariant. It is in principle possible to impose some set 
of conditions on the solutions in order to select a particular gauge for A^. However, in 
this paper we will concentrate on gauge invariant quantities and, hence, any solution 
will do. 

In the following sections we will describe how we have been able to numerically 
construct the Nahm transform of a given self-dual gauge field configuration on the torus 
by finding the solutions of Eqs. |2.15| , |2.16| using the lattice formulation of the theory. 



In the next section we will describe the numerical technique and in the following we 
will apply our construction to some explicit examples. 

3. Solving the homogeneous and inhomogenous Dirac equation 
on the lattice 



In this section we will explain how we can actually obtain the solutions of Eqs. |2.15| - |2.16 
by studying the problem on the lattice. 

All studies of the Dirac equation on the lattice have to deal with the well-known 
problem of fermion doubling. This usually means that one must either break chiral 
invariance explicitly by using Wilson fermions or maintain it at the expense of producing 
spureous lattice solutions to the equation. A middle way is represented by the use of 
staggered fermions. To put Dirac equations on the lattice we have used both the Wilson 
action, keeping the value of the Wilson parameter r free, and a naive fermion action; 
in both cases, we set the lattice fermion mass equal to zero. The second case allows for 
Weyl (2-component) spinors, while the first one requires Dirac (4-component) spinors. 
The results we report here have been obtained with Wilson fermions, so we restrict our 
discussion to this latter case. 

Our Wilson-Dirac operator reads (the superscript L stands for lattice): 

Tp L ^{n) = Ar4>(n) - \ £J (r - ^)U^n)^(n + ft) + 

+(r + 7 M )C//(n-/i)^(n-/i)] , (3.1) 

where is the gauge field, fi is a unitary lattice vector in the \i direction, and 7 M are 
Dirac matrices. To add a constant U{1) potential, as required by the construction of 
the Nahm transformation, we exponentiate the whole A fl +27rz fl I to build link variables. 
Thus, lp L z has the same form as above but with replaced by exp(— 2maz fl )U fl , where 
a is the lattice spacing. Note that the argument of the exponential is conveniently 
dimensionless, as the coordinates in the dual space have dimensions of inverse length. 

Our first goal is to solve numerically homogeneous Dirac equations in a given back- 
ground lattice gauge configuration. For this purpose, we construct the hermitian posi- 
tive operator (75 p L ) 2 , for which we expect to identify a lowest eigenspace of smooth 



modes, whose degeneracy must be equal to the one given by the index theorem. The 
elements of this space can be identified as a lattice approximation to the continuum 
zero modes. Although in general no exact fermionic zero modes will live in the lattice, 
the lowest eigenvalue of (75 p L ) 2 must approach zero with finer discretization to signal 
the recovery of continuum behaviour. All this requirements are indeed verified by our 
results, and thus we translate the problem of solving homogeneous Dirac equations into 
that of finding out these lowest eigenvectors. 

Other problem we face when computing solutions is that of keeping the right chi- 
rality for the fermion fields. Our technique to achieve this proceeds by tuning the value 
of r. There is a convenient window within the interval [0, 1] in which r is large enough 
to push up the doublers, and small enough so as to maintain approximate chirality 
for the resulting solutions. This range of values can be determined by computing with 



different r, then checking the properties of the result with different criteria |L8], |19| . A 
typical value in our computations is r = 10~ 2 . 

We are finally able to obtain vectors that are approximately chiral, the "wrong" 
half (chirality opposite to the one dictated by the sign of the topological charge of 
the gauge configuration) of the spinor being numerically negligible everywhere with 
respect to the "right" half. Above the lowest space, whose degeneracy is given by the 
index theorem, the typical structure for the lowest part of our spectra presents a first 
excited continuum level, the order of magnitude of its eigenvalue being stable with 
changing lattice volume, and a set of doubler states between this excited level and 
the zero modes level. We can recognize the doublers by the oscillating behaviour of 
their invariant densities, as well as by comparison with the known doubling symmetry 
pattern. Moreover, the evolution with r of the respective eigenvalues for the quasi-zero 
modes and the doublers is qualitatively different. The gap between the approximate 
zero modes level and the doubler levels is determined essentially by the value of r. 

To compute lowest eigenvectors of the operator (75 p L ) 2 we have used a standard 
conjugate gradient (CG) method, basically the one described in [OT (without imple- 
menting the acceleration method described in this reference), which gives eigenvectors 
of the hermitian, positive definite operator A with increasing eigenvalue by minimizing 
the Ritz functional: 

V>\?) = — — , 3.2 

< v,v > 

with v a trial CG vector. 

To check by consistency our CG method we have performed alternative calculations 
using both a Lanczos procedure and a minimization technique inspired by the cooling 
algorithm used in lattice Yang-Mills studies. The results for the zero modes in a 
number of qualitatively different gauge configurations will be the subject of future 
publication [ |T9|j . 

To solve inhomogenous Dirac equations of the form pip = <p, as required by our con- 
struction of the Nahm transformation through a local Taylor expansion, we use a vari- 
ant of the CG method, the stabilized biconjugate gradient algorithm (BiCGStab) []21j], 



whose adaptation for the purpose of inverting lattice Dirac operators is described in . 
We have used directly the implementation described in this latter reference, with minor 
operative changes. The basic idea of the construction is the following: if we apply CG 
to an inhomogenous equation of the form Av = w, for each iteration we will have a 
trial solution v^, a search direction p^ and a residue = w — Avf, the new trial vector 
is constructed as v i+ i — v $ + api, with a determined by the requirement that the new 
residue has minimum norm, and the problem is to find the new pi + \ (which fulfills 
the same condition for search directions as in "raw" CG) and r^+i with the lowest 
computational cost, which essentially means with a lowest number of operator- vector 
multiplications. The BiCG method incorporates an additional sequence of vectors 
such that r+fj = r+rj = for i = 0, 1, . . . ,j — 1. With this construction, the residuals 
satisfy 



P 3 {A)r , f 3 = PM + )ro 



where Pj are polynomials of degree < j. If the method converges, these polynomials 
can be identified as 'reductors' for the norm of the residue. 

This convergence behaviour can be enhanced as in the conjugate gradient squared 
(CGS) algorithm, which furthermore avoids the computation of the fj and makes resid- 
uals satisfy r 3 - = P?(A)r , with the same Pj. However, the convergence of CGS is 
usually far from being smooth, and the possibility exists of obtaining large peaks for 
| r^l which can plague a systematic application of the method. The BiCGStab cures this 
problem by making rj = Qj(A)Pj(A)ro, with an appropriate form for the polynomial 
Qj so as to keep a good reduction behaviour while keeping in a smoother regime the 
evolution of |rj|. Further technical details can be found in the references. 

The application of this method to our problem allows us to efficiently compute 
approximate solutions, with typical values for \<f> — far below the corrections to 
continuum quantities eventually coming from other sources, thus not adding significa- 
tive new errors to the calculations in which these vectors enter. 



4. Application to the case of non-trivial 577(2) gauge field con- 
figurations on the torus 

In order to test the applicability and stability of our method we have first studied a case 
for which an explicit analytical solution is known. This occurs for gauge configurations 
with constant field strength. The Dirac equation in this background can be solved by 
techniques which are very similar to the ones used to study small deformations around 
these gauge fields IfLOfl . An explicit example has been worked out in Ref. |p3| . For 
comparison we take precisely this case in our numerical work. 

Our gauge configuration has twist tensor defined by n 03 = n 2 \ = 2, the other 
components being zero; its topological charge is 2. The gauge potential and the field 



o 



strength tensor can be written as: 



4. 


1 

= 2 n il 


F 


= 7l li a3 



x u a 3 (4.1) 

(4-2) 

where Z M are, as defined before, the lengths of the torus, which we set equal to 1. (Note 
that this F is actually anti-self-dual). The Nahm transformation can be carried out 
analytically in this case |23|, and the result is = —A^, F^ u = F^ v = —F^ v with the 



natural gauge choice for the transformed fields. 

To compare with this result, we first put the configuration on the lattice, which 
can be done exactly by exponentiating the above form for A^ to build links. Now, 
the required numerical solutions to lattice Dirac equations can be computed, and the 
transformed field extracted from their scalar products according to our procedure. To 
be able to obtain a large set of transformed points, we have worked in a 8 4 lattice. At 
the end we get a perfect structural reproduction of the exact results, with only a slight 
error for the field value. The only nonvanishing fields for our numerical transformation 
are E 3 and B 3 , as required, and self-duality is verified almost exactly within the ma- 
chine numerical precision; when a gauge is fixed (by directly diagonalizing the fields 
at each z-point, because the gauge fixing procedure we describe below for nontrivial 
configurations does not make sense in this case), we end up with: 

E 3 = B 3 = 1.05£ 3 C , (4.3) 

with E^ the continuum field. The 5 % difference must come from the discretization 
corrections to the solutions of Dirac equations, which are sizable due to the small num- 
ber of lattice points. In any case, we can conclude that our numerical implementation 
of the Nahm transformation is working satisfactorily for this test example. 



Now let us consider another case for which no analytical solution is known. Moti- 
vated by our own interest in this type of configurations, our starting point is a self-dual 
SU (2) gauge field configuration with non-trivial twist and fractional topological charge. 
In particular, we take a solution having non-orthogonal twist given by the twist tensor: 

n 0i = (1, 0, 0) ^e ijk n jk = (1, 1, 1) (4.4) 

Given this twist, the topological charge must be a half integer value, equal to the 
action divided by 87r 2 . We are interested in the absolute minimum action solution in 
this twist sector. This has action 4n 2 and topological charge |. A configuration of 
this type, which is unique modulo space-time translations and discrete symmetries, is 
easy to obtain numerically by minimizing the lattice action in this twist sector (see, 
for example, Refs. O, [12]). The action density consists on a single lump with size of 



the order of the smallest torus length. For reasons that will be clear soon we took the 
torus sizes in the ratio l : li : l 2 : 1% =2:1:2:1. 

However, SU(N) configurations with twist are not good starting points for the 
Nahm transformation since these boundary conditions are singular for spinors trans- 
forming in the fundamental representation. There is a number of ways out of this 
problem. We will choose the following one. Instead of considering a single torus period 
we will consider more than one period in the short directions. The resulting configu- 
ration, after applying an appropriate gauge transformation, has no twist or singularity 
in the new torus of size / M = 1. Gauge invariant quantities are periodic with period 
in the 1 and 3 directions. Hence the action density consists on four lumps like those 
described previously and the total topological charge and action are 2 and 16n 2 respec- 
tively. In Fig. m we plot the action density in the x\ — 23 plane setting the other two 
coordinates to the value giving maximal action density. This is the original gauge field 
configuration to which we will apply the Nahm transformation. This will again be a 
self-dual SU{2) configuration with topological charge equal to 2. 

In the numerical construction of the Q = 1/2 lump we have worked in a lattice of 
size 12 x 6 x 12 x 6. After "glueing" the lump to itself to obtain the configuration 
explained above, the result is a 12 4 lattice, where all the subsequent computations are 
carried out. 

The first thing we did was to look at the solutions of the Dirac equation in the 
background of this configuration. Indeed on the lattice there are no exact solutions. 
What we actually did, as explained in the previous chapter, was to find the smallest 
eigenvalue of the operator (75 p L ) 2 and the corresponding eigenvectors. There are 2 
linearly independent degenerate solutions. We will choose the following procedure to 
choose a basis in this two-dimensional space (which fixes the gauge for F). Consider 
one of the points xq on the lattice where the action density peaks (there are 4 points 
of this type). Now we choose the basis of solutions of the Dirac equation such that the 
hermitian matrix M Qj g(xo) is diagonal: 

iVM* ) = (* a (*o)) + ^(* ) = (q 1 A °J • (4-5) 

It is a scalar product (sum over spin and color indices) of the 2iV-dimensional vectors 
ty a (xo). We choose Ai > A2. This condition does not fix completely the solutions, one 
is still free to multiply each element of the basis by an independent phase. We fix the 
remaining arbitrarily by imposing M a p{x' a ) to be real, where x' Q = xo + |/i is another 
of the points where the action density peaks. There is an overall choice of phase which 
will not affect our results. Our first result will be to give the density matrix M Q , ( g(x), 
which is a (gauge invariant with respect to transformations of the original gauge field) 
2x2 hermitian matrix field defined on the original torus. In Fig. |] we show the value of 
the corresponding densities M\\ and M22 in the same plane where we gave the action 
density of the configuration before. Notice that the first solution peaks precisely at 
a point Xq where the action density peaks, and was used to fix a basis in the space 
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of solutions. On the contrary the density of the other solution M 2 2 shows the same 
behaviour but displaced by a vector (0, |, 0, |) in space. Other interesting property of 
the solution is that the matrix M a p(x) turns out to be real at all space points. 

The next step was to study the solutions of the Dirac equation for other values of 
z. We use the same conditions to choose a basis in the space of solutions. The nice 
symmetry properties of the z = solutions are lost for arbitrary values of z, but the 
solution seems to show a nice smooth behaviour which makes it reliable to consider 
them a good approximation to the continuum limit ones. 

Now, we consider the solutions of Eq. p,16| , where the left hand-side is fixed by the 
solutions of Eq. 2.15 , as explained before. The solution is now unique, for a given choice 
of the direction \x and of the index a. This is in contrast with the continuum case and is 
precisely due to the fact that the lattice operator has no exact zero-modes. Fortunately, 
we are interested in F^ u which is insensitive to the continuum solution chosen, hence we 
feel satisfied with the unique solution provided by the lattice. For all values of z studied, 
the resulting ^(x, z) obtained with this procedure look pretty smooth, suggesting that 
indeed we are dealing with lattice approximants to the continuum functions. Finally 
by use of Eqs. |2.17| , |2.19| we can calculate the components of F^ u . 

The first thing to check is that the result obtained is indeed a self-dual field-strength 
tensor. For all the z values studied this was very approximately the case. In Table 1 we 
show the results for some selected values of z. The departures from self-duality were 
always found to be at most of order 10 -2 , as for the original lattice gauge configuration. 
We have actually explored points in different regions of the dual torus (periodicity in 
was checked and found to work extremely well) and always got results which indicate 
regularity SIS db function of z. 

In Fig. |3| we show a plot for the action density obtained in the plane (z2 = 1/2, z± = 
1/2), where a single peak appears. A total of 20 points in this plane have been generated 
and we used the symmetry under ~ rotations around the center of this plane to increase 
the statistics of the surface. This symmetry has indeed been checked and works to a 
high level of precision. 

The action density of the Nahm-transformed solution resembles very much the ac- 
tion density of the original configuration. When comparing one has to take into account 
that the centers (the actual maxima of the action density) of the original configuration 
do not lie exactly at a lattice point. To make a comparison which would be more pre- 
cise, we made a fit to the action density of the original lattice configuration using a few 
Fourier components in each variable and this allowed us to interpolate the information 
of the original lattice configuration to points lying in between the lattice points. In 
Fig. |] we compare the result of this fit (solid curve) to the actual data points obtained 
for the transformed field along one of the straight lines in the torus joining two cen- 
ters. We emphasize that the curve is obtained by fitting the original configuration and 
not the Nahm-transformed one. In all the points explored we always found that the 
action density of the transformed field matches quite precisely the one of the original 
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configuration. We proceeded to investigate systematically what is the exact relation 
between the original configuration and its Nahm transform by looking at the rest of 
the gauge invariant local quantities. These are the traces of the products of two spatial 
components of the electric field Tr(Ei(z) Ej(z)). We made a careful selection of the 
z points to distinguish the different possibilities. The differences between the original 
and transformed action densities in this comparison can grow up to 5-10 %, which can 
be partly due to small uncertainties in the determination of the peak centers for the 
lattice configurations, from which the direct coordinate comparison depends. However, 
we found a consistent way to recognize directions by looking at the signs of the traces 
of crossed products of fields. Finally, we arrived at the conclusion that, up to a gauge, 
the relation between the original and transformed field tensors is as follows: 

F(z) = F(x'(z)) (4.6) 
x'(z , z u z 2 , z 3 ) = (zi, -z , z 3 , -jzjj) (4.7) 

Actually, changing z — > —z is also possible since the configuration is invariant under 
this transformation. To show an explicit example of the comparison we applied the 
Nahm transformation precisely at a transformed lattice point. The result is shown in 
Table 2, where the traces Tr(Ei(x) Ej(x)) are compared for two corresponding x and z 
points. To arrive to similar results for A^, we should be able to compute and compare 
Polyakov lines. This is computationally very expensive with the present method, al- 
though some variation of it could be particularly helpful in this respect (cf. comments 
at the end of the concluding section). 



To finish this section we want to mention another application of our numerical 
methods to the study of a given gauge configuration. It comes from the study of 
fermion zero modes in the adjoint representation of the gauge group, which can be 
directly adressed with our techniques. In this case there is no need to have strictly 
periodic gauge configurations, and we can consider also twisted background fields. The 
index theorem indicates that 2NQ zero modes are to be found in this case. Furthermore, 
in contrast to the case of the fundamental representation, a new important property is 
expected: the supersymmetry present in the theory with adjoint fermions supplies 2 of 
the zero modes as supersymmetric transforms of the background (anti) self-dual gauge 
field P4] , through 

^\x) = \*t?W< > ( 4 - 8 ) 
where a is an adjoint color index in the Pauli basis, a an index in the space of zero 
modes, A, B are spin indices, u is the spinor parameter of the transformation, and 

Our lattice computations are again able to reproduce these continuum results p 



19]. In this case, the expected 2iVQ-dimensional space of lowest modes actually splits 



into NQ spaces of dimension 2, above which doubler levels appear. Each couple of 



zero modes present a geometrically distinct character, and the space related to the 
'super symmetric' zero modes can be easily identified, typically being the highest one. 
This allows to extract numerical approximations for the continuum field F^, which are 
expected to differ from the lattice values computed from plaquettes, as finite lattice 
spacing corrections enter these two quantities in completely different ways (a complete 
analysis of discretization errors in this new approach is still lacking). We have verified 
that this is indeed the case, the difference being most appreciable around absolute 
maxima and minima for the action density. However, the results for the gauge field 
coming from adjoint zero modes still compares successfully (although with larger errors) 
with the Nahm-transformed field when checking the correspondence between them. 
Further study is required in order to understand the differences in lattice corrections 
and to obtain a better way to compare the results consistently. 

5. Conclusions and future outlook 

In this paper we have presented a way in which one can numerically construct the 
Nahm transform of a self-dual Yang-Mills classical configuration discretized on the lat- 
tice. This opens the door for a more thorough study which will allow to disentangle 
the main features of this transformation. For example, one can study fixed points of 
the transformation, and relate solutions belonging to different gauge groups. This goal 
however requires an important computational effort and at the moment we are improv- 
ing our algorithm and codes to make this a feasible task. An analysis of systematic 
errors is also underway. The analytic information on the transform can be of great help 
in suggesting the appropriate choice of initial configurations and points to study. 

We have applied the Nahm transformation to a self-dual gauge field configuration 
built out of one carrying non-orthogonal twist. We found that in this case the trans- 
formation gives back the original configuration at a transformed point in space-time. 
This sort of constraint could turn out to be useful in finding the analytical solution. 

To conclude, we want to comment on a point which was left out and will be de- 
veloped in future publications. This is the fact that one can use techniques similar to 
the ones described to obtain a lattice Nahm transformation (that is, a lattice gauge 
configuration which is the Nahm transform of another lattice gauge configuration) by 
directly constructing link variables instead of continuum ones. This has a number of ad- 
vantae parallel method; for example, Polyakov loops for the Nahm-transformed 

configuration can now be easily computed from lattice links. 

The construction proceeds as follows. Let us first select a grid on the dual torus 
z^ = an^, where a is the spacing in this new lattice. The solutions ty a (x, z) of Eq. |2.15 
can now be computed in exactly the same way we have described. Once they are 
known, form the matrices: 

Vf{z) = * a+ (^ * + A) (5.1) 



where color and spin indices are summed over and jl is a unitary vector in the ^-lattice. 
These matrices transform as lattice links under a ^-dependent change of basis in the 
space of zero modes, i.e. a gauge transformation in z-space: 

V a {x, z) — ► y' a {x, z) = z)VL +Pa (z) (5.2) 

Vf — > = n ai (z)Vfn +&P (z + A) (5-3) 

but are not unitary, as required to be identified as link variables. To extract a unitary 
matrix from V M while keeping the transformation property we can decompose it as the 
product of a unitary matrix times a hermitian matrix 1 , V^z) = H ^(z)U ^(z) . Finally, a 
SU (Q) matrix has to be extracted from each U^, which in general is in U(Q). There is a 
number of ambiguities here, having to do with remaining freedoms in the 'unitarization' 
process, but these can in principle be dealt with by appealing to continuity arguments 
when going from one z point to the following, once a convenient gauge has been chosen 
(for which the method we have described in the preceding section remains valid). 

It can be seen that the naive continuum limit for the new link variables is indeed the 
continuum field in Eq. |2.1C , although a rather complicated behaviour arises for lattice 



corrections, which makes it difficult to keep them controlled. This latter point seems 
to be the main difficulty to exploit this method. Further study, and probably use of 
larger lattices is required before it can be efficiently implemented. 
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1 For Q = 2 there is an easier way: the matrix = V AI +er 2 V*c72 nas the property U^U^ = deKjA^I, 
and it is straightforward to extract a SU(Q) matrix from it. 
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Figure 1: Action density for the original lattice gauge configuration in the plane of the four 
absolute maxima. The total action has been normalized to 1. 




Figure 2: Invariant densities for the two orthogonal zero modes chosen for the original lattice 
configuration, in the plane of the four absolute maxima. Both have norm 1 in the continuum. 





Figure 3: Action density for the transformed gauge configuration in the {zq = 1/2, z% = 1/2) 
plane. The normalization is the same as in the original lattice configuration. The center in 
this plane has been displaced to the origin of coordinates for clarity, and an interpolating 
procedure has been used to obtain a smooth surface. 
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Figure 4: Comparison of action densities for the Nahm transformed field (dots) and the 
original lattice field (solid line) along a line with two absolute maxima. The line has been 
extracted from a fit to the lattice points densities. The center has been displaced to the origin 
of coordinates for clarity. 




Table 1: Values for the transformed electric and magnetic gauge fields in some selected points 
of the z-space are given. The number of significative figures is dictated by the deviations from 
self-duality appearing in the original lattice configuration. The first point corresponds to an 
absolute maximum; the second, to the maximum in a plane in which E\ = E 3 = 0; the third, 
to a point with no special character. 



z-space point 


rpa rpa rpa 


pa Da Da 


(1/2,1/4,1/2,1/4) 


(-0.049,0.012,-1.262) 
(-1.204,0.444,0.051) 
(-0.436,-1.182,0.005) 


(-0.049,0.012,-1.262) 
(-1.212,0.447,0.051) 
(-0.436,-1,182,0.005) 


(1/2,1/2,1/2,1/2) 


(0.000,0.000,0.000) 
(0.218,-0.705,0.000) 
(0.000,0.000,0.000) 


(0.000,0.000,0.000) 
(0.225,-0.725,0.000) 
(0.000,0.000,0.000) 


(3/4,1/2,1/2,2/3) 


(0.031,-0.001,-0.262) 
(-0.311,-0.166,-0.038) 
(0.076,-0.287,0.010) 


(0.032,-0.002,-0.262) 
(-0.316,-0.169,-0.037) 
(0.076,-0.287,0.010) 



Table 2: Values for the traces of products of electric field components for corresponding 
x and z points: x = (0.312,0.584,0.596,0.766), z = (0.584,0.688,0.766,0.404). x cor- 
responds to a lattice point, selected such that clear hierarchies among the invariants are 
established. 





Xq 


Zq 


Tr(E!(x) E^x)) 


0.418 


0.460 


Tr{E 2 (x) E 2 (x)) 


0.465 


0.488 


Tr(E 3 (x) E 3 (x)) 


0.383 


0.412 


TriE^x) E 2 (x)) 


0.033 


0.031 


Tr{Ei{x) E z (x)) 


-0.027 


-0.027 


Tr(E 2 (x) E 3 (x)) 


0.024 


0.025 



on 



